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Abstract 

The propagation of light-pulse with negative group- velocity in a nonlinear medium 
is studied theoretically. We show that the necessary conditions for these effects to be 
observable are realized in a three-level A-system interacting with a linearly polarized laser 
beam in the presence of a static magnetic field. In low power regime, when all other non- 
linear processes are negligible, the light-induced Zeeman coherence cancels the resonant 
absorption of the medium almost completely, but preserves the dispersion anomalous and 
very high. As a result, a superluminal light pulse propagation can be observed in the sense 
that the peak of the transmitted pulse exits the medium before the peak of the incident 
pulse enters. There is no violation of causality and energy conservation. Moreover, the 
superluminal effects are prominently manifested in the reshaping of pulse, which is caused 
by the intensity-dependent pulse velocity. Unlike the shock wave formation in a nonlinear 
medium with normal dispersion, here, the self-steepening of the pulse trailing edge takes 
place due to the fact that the more intense parts of the pulse travel slower. The pre- 
dicted effect can be easily observed in the well known schemes employed for studying of 
nonlinear magneto-optical rotation. The upper bound of sample length is found from the 
criterion that the pulse self-steepening and group-advance time are observable without 
pulse distortion caused by the group-velocity dispersion. 

PACS numbers: 32.80.Bx, 42.50.Gy, 42.50.Gt, and 42.25.Bs 

1 Introduction 

The propagation of light pulses through the media with a normal dispersion is adequately 
described by the group velocity, which in this case turns out to be, in fact, the velocity of 
signal transfer. However, the group velocity loses its role of signal velocity in the region of 
anomalous dispersion, where it not only exceeds the vacuum speed of light, but at certain 
frequencies becomes negative, whereas, even when the group velocity is superluminal, no real 
signal can be transferred faster than the vacuum velocity of light Nevertheless, contrary to 
the common assertion p, |^ that superluminal group velocity is not a useful concept, Garrett and 
McCumber 0] have shown that in the case of superluminal propagation of sufficiently smooth 
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pulses, such as Gaussian wave packets,the motion of the pulse peak is correctly described by 
the classical expression of group velocity. The effect of superluminality is that the emerging 
pulse has essentially the same shape and width as that of the incident wave packet, but its 
peak exits the cell before the incident pulse even enters. This process can be understood in 
terms of superposition and interference of traveling plane waves that have a distribution of 
frequencies and add up to form a narrow-band light pulse. In a dispersive medium the speed of 
each plane wave is determined by the refractive index at that frequency. Around the resonance 
in the region of anomalous dispersion the refractive index varies so steeply and different plane 
waves interfere in a way that the point of maximal constructive interference is localized at the 
exit face of the sample. The plane waves appear in the medium much earlier and, hence, the 
exit pulse is formed long before than the peak of incident pulse enters the medium. When the 
group velocity is negative the exit pulse splits into two pulses, one of which moves forward in 
the vacuum, while the second one propagates backward in the medium and is canceled by the 
incoming pulse at the entrance of the cell (see Sec. 3). The super luminal propagation can be 
phenomenologically represented as reshaping of the pulse, most of which is absorbed, leaving 
only a small pulse in the leading edge. Eventually, the outgoing pulse moves superluminally 
and saves the shape of initial pulse. This picture breaks down, if the envelope of incident wave 
packet has singularities, for example, step-type modulations, that always moves with the light 
velocity c. All these results demonstrate that for smooth pulses with continuous derivatives the 
group velocity is physically meaningful even when it is negative. But it cannot be interpreted 
as a velocity of real signal, which could propagate no faster than light These predictions 
have been confirmed experimentally using various schemes, including the anomalous dispersion 
near an absorption line P, |^, nonlinear gain [§, active plasma medium 0] and the tunneling 
of one-photon wave packet through a barrier [|1^ . 



An interesting model for a transparent anomalous dispersive medium and a light pulse 
propagation with negative group velocity has been proposed by Chiao and co-workers |TT], |T^. 
The model is based on the destructive interference between the two close-lying gain transitions. 
In this scheme L.Wang and co-workers have recently observed a pulse-advance shift of 62ns in 



cesium vapor [T^ 



The results of papers [^-pf, 1131" based on linear response of the medium where 



the negative group velocity always exists within the absorption lines M and outside gain lines 



TT|]-0]. The similar results have been obtained, when a weak probe beam propagates with 



an effective index of refraction induced by a nonlinear interaction of the medium with a strong 
pump field. Under these conditions both the probe light ultraslow propagation Il5|-[l^ at the 
electromagnetically induced transparency (EIT) and the negative "group-delay" propaga- 
tion associated with electromagnetically induced opacity [jl9|, ^ have been recently observed 
experimentally. However, from the point of view of the dependence of optical properties of 
medium on probe field intensity these systems may be clearly considered as effectively linear 
optical media with an extraordinary large dispersion. 

Intriguing questions arise whether the nonlinear superluminality occurs in a medium, the 
refractive index of which is modified by the propagating light pulse. How this effect is pro- 
nounced, and what is the form of superluminal pulse emerging from such a nonlinear medium? 
It seemed, that the answer to these questions could be easily obtained in a two-level atom, 
which is a basic model for the study of superluminal propagation in hnear media p|-p[],p!0| . 



Notice that this system has been also employed for investigation of nonlinear slow optical pulse 
propagation in media with normal dispersion |2T|, |22[. However, a two-state system is unfit 



for observation of nonlinear superluminal effects. The problem is to observe the pulse super- 
luminal propagation in a medium with an intensity-dependent refractive index, but without 
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accompanying nonlinear processes, which lead to pulse reshaping as well. Since a two-level 
atom displays anomalous dispersion within the absorption line and the nonlinearity in the dis- 
persion is achieved at the intensities that are much larger than the saturation intensity (see 
Appendix), nonlinear effects such as self-focusing (or defocusing), self-phase modulation etc., 
cannot be neglected. Meanwhile, the systems, that are almost transparent and at the same 
time show large nonlinear dispersion, have been basic models for theoretical and experimental 
studies of resonant nonlinear magneto-optical effects caused by a laser-induced coherence (see 
the recent review by Budker et al. |2^). The simplest system displaying these properties is a 



three-level A-atom with ground-state Zeeman sublevels, that interacts with the right (RCP) 
and left (LCP) circularly polarized components of optical light-pulse in the presence of weak 
magnetic field, which is applied along the direction of light pulse propagation (Fig.^). Due to 
the nearly maximal Zeeman coherence, large nonlinear susceptibilities are created at the reso- 
nance frequency in analogous to those observed in EIT experiments (for relationship between 
EIT and nonlinear magneto-optics see fl^ ) . This leads to very narrow magnetic resonances p4 1 
or enhances significantly the nonlinear Faraday signal [^. However, the dynamics of resonant 
light propagation in coherent nonlinear media has not been studied yet. This analysis is carried 
out for the first time in the present paper. For the Faraday configuration in Fig.|I| we show that 
the light-induced Zeeman coherence preserves the anomalous dispersion of the medium, while 
it suppresses the absorption at the resonance line center. This results in a negative nonlinear 
group- velocity at the values of light intensity, that are many times smaller than the saturation 
intensity of the corresponding atomic transition. In this scheme the superluminal effects are 
revealed in two ways. Firstly, similar to linear media, the peak of outgoing pulse exits the cell 
before the peak of incident pulse with a narrow spectral width enters the medium. Secondly, 
the pulse velocity is nonlinear in light intensity, the more intense parts of the pulse travel 
slower. As a result, the self-steepening of the trailing edge of pulse takes place in contrast to 
the shock-wave formation in media with normal dispersion I^T], In fact, the steepening 
effect is the most striking manifestation of anomalous dispersion, if one takes into account that 
usually the residual absorption causes a misconstruction of the pulse advance. 

The paper is organized as follows. In the next section we derive equations for the evolution 
of the atomic density matrix and for the field transmitted through the medium. Basing on 
the analytical solution for the density matrix we give the results of numerical calculations for 
the pulse propagation with nonlinear negative group velocity in Sec. 3. Our conclusions are 
summarized in Sec. 4. Finally, in the Appendix we calculate the nonlinear group-velocity of 
light pulse propagating in the Doppler broadened medium for the cases of two-level atom and 
three-level A-system. 



2 Nonlinear negative group velocity in a medium with 
Zeeman coherence 

It is well known that in linear media the anomalous dispersion always occurs within an ab- 
sorption line ||l|] and leads to negative group velocity for the peak motion of sufficiently smooth 
pulses [Q. However, as can be seen from Fig.^ for the case of a two-level atom, in usual mate- 
rials the light pulses experience very large absorption in the vicinity of sharp atomic resonance 



that prevents a doubtless observation of high anomalous dispersion ||2^ . On the other hand, 
the light-induced Zeeman coherence in a three-level A-atom changes essentially this absorption- 
dispersion relationship, making the medium almost transparent, while keeping its dispersion 
still very large (Fig.^. Moreover, as it is apparent from Fig.^, there is a frequency range, where 
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the dispersion on both transitions |1) — > |3) and |2) — > |3) is the same and anomalous. This 
means that the RCP and LCP components of hnearly polarized laser field travel with the same 
negative group velocity. As the latter depends on the hght intensity, the pulse advance, being 
the minimum at the pulse peak (see below), is different for different parts of the pulse compared 
to the same pulse traveling the same distance in a vacuum. As a result, the trailing edge of 
pulse should show self-steepening that can be tested in an interferometric experiment similar 



to that performed by L. Wang and co-workers 13 



2.1 Basic equations 

The atomic system depicted in Fig.[l| interacts with the fields of RCP and LCP components £"1^2 
of linearly polarized laser light with the frequency uj and wavenumber k propagating along the 
z axis. In the presence of a static longitudinal magnetic field B these components,which are 
working on the transitions |1) |3) and |2) — > |3), respectively, are detuned from the atomic 
resonances by 

Ai,2 = w - ^31,2 =F A w = £ifA, (1) 

c 

where Uik is the frequency difference between the levels i and /c, and A = gfiBB/h is the 
Zeeman shift, which is induced by the magnetic field for the ground state sublevels |1) and |2) 
having magnetic quantum numbers =Fl. Here g is the ground-state gyromagnetic factor and hb 
is Bohr's magneton. In Eq.([^) we have included the Doppler shift -v for atoms moving with 
frequency Vz = v. 

The interaction of fields Ei^2 with a single atom is described by the Hamiltonian: 

Hint = hAicrn + /iA20-22 - ^(^^10-31 + ^^20-32 + h. c). (2) 

Here aij = \ are atomic operators. The Rabi frequencies of field polarization components 
are defined as Qi^2 = /^-£'i,2/^; where the dipole moments of the transitions |1) —>■ |3) and 
1 2) — > 1 3) are taken equal: = ^23 = Z^- 

The time evolution of the system's density matrix p obeys the master equation 

= -^[^mt ,P] + Ap, (3) 

where the matrix Ap accounts for all atomic relaxations. In further, we neglect the collisional 
broadening of optical transitions compared to the spontaneous decay rate, but we take it into 
account for the ground-state coherence. This approximation is valid for the gas pressure below 
lO^^cm"^ at room temperature that is assumed hereafter. Thus, in our model the longitudinal 
and transverse optical relaxations of atomic Bloch vectors are determined by the process of 
spontaneous emission from the upper level |3) to the lower levels |1) and |2) with equal rates 
7i = 72 = 7; and outside of the three-level A-system with the rate 70. The Zeeman coherence 
decay rate 7c is determined by the collisional and atomic time-of-fiight broadening. Thus, the 
phenomenological relaxation matrix Ap has the form 

/ 7P33 -7cPl2 -Ppis \ 

Ap= -7cP21 7P33 -rp23 , (4) 

V -Fpsi -Fpsa -(7o + 27)p33/ 
where 2r = 70 -|- 27 is the total spontaneous decay rate of the level |3). 
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The evolution of the slowly varying field amplitudes E12 along the z axis is determined by 
Maxwell's equations 



(5) 



where Pi^2{z,t) are the field induced polarization, real and imaginary parts of which are re- 
sponsible for the dispersive and absorptive properties of medium at the transitions |1) |3) 
and 1 2) |3), respectively. The polarizations are calculated from the Fourier transform 



P^{z,t) 



1,2, 



(6) 



where u = uj — u is the deviation from the field carrier frequency uj, and Pi(z/, z) are expressed 
by the corresponding susceptibilities 



as 



(7) 
(8) 



= Xi{^, z)Ei{iy,z), 

with being the atomic number density. 

We consider the case of smooth incident pulses with a duration much larger than all relax- 
ation times of the medium. Then, the limited bandwidth of field Fourier amplitudes Ei{v,z) 
allows us to approximate Xi{''^i ^) by the first few terms of the Taylor series: 



X{^, z) = x{.^,z) + 



dx 
du 



v=0 



2 5z/2 



v=0 



The substitution of Eq.@ into Eq.(|) yields 



Pi(z, t) = Xii^, z) ■ Ei{z, t) + i 



■ dXi 



dui 



dE, 1 



dt 2 duj"^ 



+ 



Then, from Eq.(^) we have 
\ dz Vi dt 



Ei{z,t) = 2m 



UJ 



Xi{uJ,z)Ei - 



2 dco^ 



(9) 



(10) 



(11) 



where group velocities Vi {i = 1, 2) for the two polarization components are introduced in the 
form 



OJ 



1,2. 



(12) 



This expression coincides with the classical formula of group velocity Vgr = c 3?e 



n + u 



-1 



if one takes into account that the complex index of refraction of the medium n = 1 + 2'kx and 
3f?ex << 1- The second term in rhs of Eq.(|l^) is responsible for the pulse distortion due to the 
group-velocity dispersion. In the next section we will show that this term is negligibly small 
for the length of the atomic sample L, which is restricted by the requirement of smoothness of 
the pulse. 

Since in this paper we restrict our attention to the superluminal pulse propagation, we need 
only the equation for intensities Ii{z,t) = ^\Ei{z,t)\'^. From Eq.(|Tl|) we have 



9i + ^9^r^("'') = -"^^^^"'')' 



(13) 
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where the nonhnear absorption coefficients at the frequency a; are defined as 

ai = bmxi = ^^m-— . (14) 

Here £o = ^ind a = 47rco'/i^(/;,cr)~-'^ is the resonant absorption cross-section. 

The advance time Tad of the field Ei{z,t) propagating through a medium of length L com- 
pared to the pulse passing the same distance in a vacuum is given by 



Vi J c ouo 



(15) 



Thus, at the frequencies where the derivative |^ takes a large negative value the pulse advance 
is significant. 



2.2 Solutions 



We solve Eq.(|T3|) assuming that the incident laser beam is a Gaussian pulse 

I{z,t) = Iin{t-z/c) = Ioexp[-{t-z/c)yT% z<0, (16) 

with a temporal duration r ^ T^^. 

During the propagation of pulse through the medium the group- velocity dispersion and 
self-steepening change the pulse width that violates the approximation of truncated Taylor 
expansion Eq.(|IO|). Our aim is to observe the pulse nonlinear steepening under superluminal 
propagation without pulse spreading. The source of the latter is the derivative d'^x/d^^'^, which 
enters in the second term in the rhs of Eq. ([TI|) . Since this term is also proportional to the 
second time-derivative of the field d'^E/dt'^, the contribution from the group- velocity dispersion 
dramatically increases at the distances, where a steep front is formed due to the dependence of 
group velocity on light intensity. In order to avoid this complication, as well as to use rightly 
the steady-state solution of Eqs.(0), we restrict the upper bound of the sample length by a 
value for which the pulse broadening caused by the self-steepening does not exceed the medium 
absorption line-width F, i.e. the condition 

dl{z,t) 



dt 



<TI{z,t), 0<z<L, 



(17) 



should hold everywhere inside the medium for any point on the pulse envelope. Then, by the 
condition ([l^), the spreading of the wave packet caused by the group- velocity dispersion can 
be neglected up to the distances restricted by the value 



2 



18) 



Eq.(|T8D follows from the smallness of the second term in the rhs of Eg . (|TT|) . Here we have 
replaced the derivative d'^x/duj'^ by its maximal value taking into account that x is an implicit 
function of z through its dependence on pulse intensity. 

The condition (jl^) allows us to use the steady-state solution for atomic density matrix 
elements. This solution is obtained from Eq.(^ in the form 

P3i 



3?e 



P31 



qrA'^6p/D, i=l,2, 
-A{qeA + n^)Sp/D, 
-A{qeA-n^)6p/D, 



(19) 
(20) 
(21) 
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where Q"^ = jf^ip = |fi2p- It is derived under the assumptions of A, 7c -C F and > Qq ^ 
T'jc, where Qq is the Rabi frequency corresponding to the peak value of pulse intensity. The 
population difference between the ground-state Zeeman sublevels and the upper state |3) is 
5p, which is 5p ~ I for equal decay rates 7o = 7 and a weak magnetic field. In the linear 
limit < Ar) these expressions coincide with that for two-level atom presented graphically 
m Fig.H In the inverse case of ^ AF the medium becomes almost transparent near the 
resonance e = 0, while its dispersion remains highly anomalous (Fig.^). This result differs 
from the well-known fact established for linear media that the optical dispersion of transparent 
medium is always normal |2^. The latter directly follows from Kramers- Kronig relations, 
that assume linearity and causality in the response of the medium. In our case this response 
is nonlinear and the classical Kramers-Kronig relations are no longer valid. The numerical 
solution of Eq.(|]) in steady-state regime is shown in Fig.||. Our results given by Eqs. (p!9D -(p2D 
are in excellent agreement with these calculations. 

In further, we need to take into account the Doppler broadening by averaging the matrix 
elements pij over the atomic velocity distribution. Using the results of calculations performed 
in Appendix we finally obtain the group velocity Vgr and the absorption coefficient a of two 
polarization components of the field in the form 



1 1 

" - £o' Al, 'b{b + iy 



where b is given by Eq.(^).As it follows from Eq.(|3D, the group velocity index Ug reaches its 
maximum at the peak value of light intensity. As a result, more intense parts of the pulse move 
slower that leads to the steepening of the trailing edge of the pulse. Also, notice that in the 
absence of external magnetic field, when A = 0, the susceptibilities Xi vanish due to complete 
population trapping and the pulse moves with the light vacuum velocity c. 

The results of numerical integration of Eq.(^) for pulse propagation using Eqs. (PHj) and 
( p4D are presented in the next section. 

Now, having the solution for susceptibilities Eqs. (|A.6|) - (|A.8|) , we can find the upper bound 



of the cell length from the conditions (|1^) and (|T^). The maximal value \d'^x/duj^\^ax ~ 
N fi'^6p/4hA%r occurs at fi^ = AF. Substituting it into Eq.(|l|) we have 

where /x^ is expressed by the radiative decay rate 7 = 32tt^p'^ /3hX^ and A is the wavelength of 
atomic transition. 

The condition (|TB]) can be rewritten in the form 

Ar « i, (26) 

where 

Ar = 2nujT 



-. m 
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is the amount of spreading experienced by a pulse, when travehng through a medium with the 
length L. Thus, condition ( p6D ensures the truly observable steepening of the backward edge 
of light-pulse without its distortion. 

In order to find the upper bound for the cell length from (|17|) we consider the situation 
where the absorption is absent. In this case the reduced equation ( |l3l) has simple wave solution 
in the form 

I{z,t) =hn{t- Z/Vg), (28) 

where Iin{z,t) is given in Eq.(|T6|). Since at the trailing edge dl/dt < 0, from Eq.(|27|) and (|T^ 
follows 

^ - - TTp ^^^^ 

Note that Eq. (|29| ) offers the most stringent restriction at the peak intensity where b is close 
to the unity. 

The last criterion is that At should be small as compared to the group-advance time Tad, 
At <^ Tad- Using Eqs.(|I3|) and (^) and solutions Eqs.(^3|-^) we obtain 

m^6p 26 + 1 



AAl b^{b + l) 



< 1 (30) 



One can easily show that this condition fails only for sufficiently small Vl ~ V AT, so that 
there is no reason to consider this complication here. 

It must be also noted that for the group-advance to be observable, it is sufficient that Tad 
is large with respect to the resolution time of interferometric detection, being at the same time 
much smaller than the pulse width r~^, as it was in the recent experiment [O. 



3 Numerical results and discussion 

Here we give results of numerical calculations based on the equations derived in the preceding 
section and discuss the relevant physics. 

In our simulations for the model of atomic system we have chosen parameters corresponding 
to the Di absorption line of ^^Rb. It is worth noting that a four-state system has been employed 
in |]25[ for investigation of nonlinear magneto-optical effects in rubidium.lt has been shown that 



this simplified model, with well-chosen parameters, represents the qualitative physics quite 



well [25|. In our case the total decay rate of the upper level |3) is chosen 2r = 27r x 5 ■ 10 s~ , 
the wavelength of optical transitions A = 800 nm, and the inhomogeneous Doppler broadening 
width A/) — 50 r. For the atomic number density = 2-10^^ cm'^ we have io — 1-25 ■ 10~^ cm. 
The external magnetic field is taken B = 30 mG, which corresponds to A = 0.002 F for the 
ground-state level of rubidium. The peak value of Rabi frequency of the field is chosen to obey 
the condition fll = SAA^, i.e. Qq = 0.6 F (/q = 5 mW / cm"^) , which is below the saturating 
value for a two- level atom Qs = 2F. The Zeeman coherence relaxation rate 7c is attributed 
to the coUisional (in the presence of buffer gas) and time-of-flight broadening and is taken 
7c < 10~^F. Then, for initial pulse duration r = and population difference Ap ~ 0.25 the 
upper bounds of the sample length are estimated from Eqs.(p5D and (^) as Li ~ 4 ■ 10"^ cm and 
L2 — 400 cm, respectively. This means that in the cell with a length of a few centimeters the 
superluminal effects including the pulse self-steepening can be easily observed with negligible 
pulse distortion. 

In Fig. 4 we plot the absorption of the pulse and the group-advance time calculated from 
Eqs.(0) and (plSj) as a function of the pulse intensity for the set of parameters given above. 
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It is apparent that the pulse peak is advanced with respect to the same pulse traveling in the 
vacuum approximately fifty times, that for L = 6 cm corresponds to Tad — 10 ns. The advance 
time increases towards to the pulse wings, and it is about 70 ns for the intensity / ~ 0.1 Jq. 
Since the absorption is strongly reduced around the pulse peak, where the EIT conditions 

r]^> AAB,7c(A^ + r), (31) 

are established, the group-advance time should be easily observed by nanosecond resolution 
time. 

As an illustration, in Fig.^ we give a space-time plot of solution to Eq.(|T3|) for the initial 
Gaussian pulse with duration r ~ IF"^. The group-velocity is negative in the medium. This 
effect appears as a "back-on-time" motion of the pulse. The formation of sharp back front of 
the pulse to the end of the cell is well visible (Fig]^). 

In order to reconstruct the superluminal propagation in details we present in Fig.^ the 
pulse motion in the form of time-sequence graphs. At first, before the incident pulse enters 
the medium, a weak pulse is created at the output. Then, this peak increases in time and 
eventually splits into forward- and backward-moving pulses with very sharp trailing edges. The 
first of them continues to move in a vacuum with velocity c, and after a long time it represents 
the final pulse. The second pulse exhibits, however, an unusual behavior. It increases while 
traveling through the absorbing medium. Nevertheless, here there is no inconsistency with the 
energy conservation law. Firstly, the velocity of energy transport defined as ve = S/W, where 
5* and W are the energy flux and fleld energy densities averaged over a period, is always less 
than c, as is shown in Secondly, the increment of the medium energy in the presence 



of light field is proportional to d{ujn)/duj ■ \E\ and is negative for anomalous-dispersive 



media. Obviously this energy cannot be stored in the medium and it is transformed into the 
backward-moving wave. However, this takes place until the peak of incident pulse enters the 
medium. After that, the real absorption of the light pulse results in a cancellation of backward 
wave at the entrance of the cell. 

In Eq.(pr3[) the nonlinear absorption and intensity-dependent group- velocity are responsible 
for the reshaping of pulse. In order to demonstrate the infiuence of the two mechanisms 
separately we solved Eq.(|l^) in a hypothetical case of nonlinear absorption, assuming that the 
group- velocity is independent on the intensity of light. Setting b=0 in Eq . (|23|) , we calculated a 
propagation of Gaussian pulse with the same set of parameters. The results depicted in Fig.|^ 
show that in this case the final pulse is built up much earlier than in the case of nonlinear 
group- velocity and, what is most important, it preserves the shape and width of incident pulse. 
A slight decrease of absorption is observed also. 



4 Summary 

In this paper we analyzed the superluminal propagation of light-pulse in a nonlinear medium 
and showed that necessary conditions for these effects to be observable are realized in a three- 
level A-system interacting with a linearly polarized laser beam in the presence of a static 
magnetic field. It is highly important that the nonlinearity of refractive index of the medium 
arises in low power regime, when all other nonlinear processes are negligible. Our results 
demonstrate that the light-induced Zeeman coherence cancels the resonant absorption of the 
medium almost completely, but it preserves the dispersion anomalous and very high. As a result, 
the transparent propagation of light pulse with a negative group-velocity, which is nonlinear in 
the light intensity, is predicted. This leads to the formation of extremely sharp trailing edge 
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and to the lengthening of the leading edge of the pulse. Such a behavior is inverse to the case 
of nonlinear propagation of the pulse in a medium with normal dispersion. The predicted effect 
is the most striking manifestation of superluminality and it can be easily observed in the well 
known schemes that have been used for studying the nonlinear magneto-optical rotation. 

The discussion of many questions remained behind the scope of this paper. In particular, 
we have not given due attention to the proof of the fact that the predicted effect does not 
violate the causality. Here we note only that the proof is carried out in the same manner as in 
previous studies of superluminal propagation ||^. Further, for real atomic systems the complete 
energy state description of the multilevel structure should be included. This question requires 
a careful examination and it is in current study, the results of which will be presented in the 
next publication. 
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Appendix 

Here we calculate the nonlinear group-velocity of light pulse propagating through a Doppler 
broadened atomic medium. The cases of two-level atom and three-state A-system are consid- 
ered. In order to obtain the simple analytical expressions, when averaging the group velocities 
over the atomic thermal motion, we approximate the usual Gaussian distribution of atomic 
velocity with Lorentzian function W{v) = A£i/7r[A|, + {kv)"^], where is full width half 
maximum of Doppler broadening. 

In the case of two-level atom the imaginary and real parts of the susceptibility of medium 
are given by 

3m X = Nij^T/h{T^ + + 2n^) , (A.l) 
gfJex = Nfi^e/hiv^ + + 2n^) , (A.2) 



where e = Uat — u + kv, and Q is the Rabi frequency of the field, which couples the upper and 
lower levels of the atom. Straightforward calculation of the velocity-averaged susceptibilities in 
the first order of Au = uo — ujat gives 

nbo{bo + Ad) 

i^ex) = .. Au, (A.4) 

n{bo + AdY 

where 6o = (r^ + 2^2^)^. For the group- velocity of light pulse from Eq. (|A.4|) we find 



(A.6) 



Similarly, for the case of three-level A-atom, from Eqs.(p!9D-(p2D we have 

26+1 



AAlh^{h+l 



-Au' + ... 
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^ hAl 6+1 



Atu^ + ... 



qbA 6 + 1 2gAA|5 62(6+1)2 



, (A.7) 



(3fJex2) = ^exi{n'^-n'), (A.8) 

where 

AAd 

Correspondingly, for the group velocity of polarization components of the field is obtained 



in the form of Eq.(^). The comparison of Eqs.( |A.5D and (|23|) shows that, for the case of 
two-level atom, the dependence of Vgr on light intensity appears at > A/? ^ F, whereas in a 
three-level A-system it is achieved at much smaller values of > (AA^))^. 
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Figure 1: Open A-configuration for studying the nonlinear superluminal light propagation. a± 
are components of a linearly polarized laser field with Rabi frequency Vl coupling the level |3) 
with |1) and |2), respectively. The Zeeman shift of ground-state sublevels Ai 2 = =fA is induced 
by the longitudinal magnetic field B. Radiative decay from |3) to |1) and |2) goes at rate 7 
and the outside with the rate 70. 
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Figure 2: Real {x') and imaginary {x") parts of the susceptibility of two- level atom as a function 
of the atom- field detaining e — uj — cuat in units of F. 
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Figure 3: Imaginary {x") and real {x[) parts of susceptibilities at a+ and (7_ transitions in a 
A-atom as a function of £ = co — LOat for = 0.3r, A = O.Oir, 7c = lO^^F. Doppler broadening 
is neglected. 




Figure 4: Absorption a (in units of £o^) group-advance time Tat (in units of L/c) as a 
function of light intensity. The parameters are given in the text. 
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22 

Figure 5: Space-time image of Gaussian pulse motion with initial pulse duration r — ir~^. 
Times arc given in units of and distances in cT'^. The medium extends from 2; = 10 to 
16. The peak of the pulse is at 2; = —7 at t = 0. The group velocity and absorption coefficient 
are calculated for the parameters given in the text. 
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Figure 6: Diagrams of Gaussian pulse propagation corresponding to sequential instants in Fig. 5. 
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Figure 7: Same configurations as for Fig. 6, but in the hypothetical case of nonlinear absorption 
without intensity dependence of group- velocity. 
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